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We report a solution of the one-dimensional Schrodinger equation with a hyperbolic double-well 
confining potential via a transformation to the so-called confluent Heun equation. We discuss the 
requirements on the parameters of the system in which a reduction to Heun polynomials is possible, 
representing the wavefunctions of bound states. 
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I. INTRODUCTION 

A small number of exact solutions to the Schrodinger equation were obtained historically in the genesis of quantum 
mechanics. 1 5 More recently, other exactly-soluble systems have been found by both traditional means 6,7 and via 
the factorization techniques of supersymmetric quantum mechanics (SUSY)££ Conventionally, analytical solutions to 
the Schrodinger equation were found via a reduction to a hypergeometric equation, 10 an equation with three regular 
singular points, however recently a new solution via Heun's differential equation^"— has been reported^ 

Heun's equation, a Fuchsian equation with four regular singular points, was initially studied by the German math- 
ematician Karl Heun in the late 19th century^ It has several special cases of great importance in mathematical 
physics, namely the Lame, Mathieu and spheriodal differential equations^ However, it is only recently that its use in 
physics has become increasingly widespread, 18 with its solutions being used in works ranging from quantum rings 19 
to black holes 3 

In this work, we report a class of confining potentials that can be transformed to the confluent Heun equation 
from the one-dimensional Schrodinger equation. One case, of a hyperbolic double-well, allows one to reduce the 
solution to simple polynomials for special values of the system parameters. The double-well problem has been studied 
extensively^"— and is of continued interest due to its importance as a toy model, from heterostructure physics^ to 
the trapping of Bose-Einstein condensates. 25 

The rest of this work is as follows. We give a full derivation of our solution in terms of confluent Heun functions 
in Sec. HH along with details of the power series solution. In Sec. Hill we provide the eigenvalues of our hyperbolic 
double- well problem, along with examples of the wavefunctions found for the first few states. We draw conclusions in 
Sec. El 

II. REDUCTION TO A CONFLUENT HEUN EQUATION 

The stationary one-dimensional Schrodinger equation for a particle of mass m and energy E in a potential V(x) is 

We consider the following class of confining potentials, shaped by the depth and width parameters Vq and d 



smh p (x/d) 
cosh 6 (x/d) 



V{x) = -V =£!% , p = 2,4,6 (2) 



such that with p = 4, which we display in Fig. [TJ one finds upon substitution of V(x) 

d 2 , , x / ,o T9 sinh 4 (z)\ , . . 
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where e = 2mE/h 2 , Uq = 2mVo/h 2 and z = x/d. Upon making the change of variable £ = 1/ cosh 2 (z), such that the 
domain — oo < x < oo maps to < £ < 1, we find 

£ 2 (! - o|^(o + £(i - |o^(o + \ K + w 2 £(i - o 2 ) m = o- (4) 
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FIG. 1: Plot of the hyperbolic double- well potential under consideration, Eq. J2} with p — 4. 



Noting that at large £ Eq. dU reduces to 



d 2 



U d 2 



2 m = ^m, 



with solutions in the form 



leads us to choose the ansatz solution 



HO = e ±f <S a = -dy/Uo, 
m = e f « /(C), 



(5) 
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which yields upon substitution into Eq. (jl]) 
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where /3 = —id^/e. Undertaking a peeling-off procedure with /(£) = <^/ 2 ?/(£), we find 



w(0 = o, 



where 



a = -dy/U^, P = -idy/e, 7 = /i = ^ + 2) + 2a/3 - /3(/3 + 1)) , v = ^ (a + /3(/3 + 1)) . 

This is Heun's confluent differential equation.— It has as a solution around the regular singular point £ = given by 
the confluent Heun function 14 



Hc(ot,P,K,S,rj,£) = ^ v n (a, /?, 7, 5, 77, 0f\ radius of convergence |f| < 1, 



(10) 



n=0 



where 



5 = LL + v-^{p + 1 + 2) = hj d\ v = ^(p + 1) - n~±(p + 1 + fr) = \(1 - (E + U )d 2 ), 

and the coefficients v n are given by the three-term recurrence relation 

A n v n = B n v n -i + C n v n _2, with initial conditions v-i =0, = 1? 



(ii) 
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where 

A n = l + ^, (12a) 
n 

B n = l + i(/3 + 7 -a-l) + ^|r/-i(/3 + 7-«)-^ + ^}, (12b) 

C B = 4f*+to +n _iV (l 2c) 
n z \a 2 J 

Thus we have found the following solution to the Schrodinger equation Eq. (j4]) 

m = ^ /2 H c (a, f3, 7 , <J, »?, 0- (13) 

There is a general theorem of quantum mechanics that tells us when the Hamiltonian commutes with the parity 
operator, we should expect to find solutions of the Schrodinger equation that are either symmetric or antisymmetric 
(at least for non-degenerate states). We can see from Eq. ([T3]) that our solutions are wholly symmetric, as one would 
anticipate from the initial variable change to £ = 1/ cosh 2 (x/<i). Therefore, we seek an antisymmetric solution using 
the more convenient odd variable ( = tanh(#/d)j2£ and we find 

HOa = C(l - C 2 f /2 e~^ 2 H c (-a, - 7) /?, S, V + a 2 /4, C 2 ). (14) 

In the next section, we shall look at the instances in which one can reduce the solution to a polynomial in £ or £. 
Please note, a derivation similar to the one above can be carried for single well case p = and also in the second 
double- well case of p = 2. However, in these two cases the confluent Heun function cannot be reduced to a polynomial. 



III. BOUND STATES IN A HYPERBOLIC DOUBLE- WELL POTENTIAL 



To reduce a confluent Heun function to a Heun polynomial of degree N we need two successive terms in the 
three-term recurrence relation Eq. (fTTj) to vanish, halting the infinite series Eq. ([TO]) . This requirement results in two 
termination conditions, which both need to be satisfied simultaneously 14 

- + ^±l + N +1 = 0, (15a) 
a 2 

A iv+1 (/i) = 0. (15b) 

The first condition arises from ensuring CW+2 = 0. The second condition, which can be written as a tridiagonal 
determinant (please see Appendix [A]) , arises from ensuring vn+i = 0, such that it follows from Eq. ([TT]) that all 
further terms in the series vanish identically. In our case, the first termination condition Eq. (|15a|) allows us to find 
the following eigenvalue spectra 

£ n = —£p ( 3 + m ~ dVu^) 2 , U d 2 > (3 + 4A0 2 , (16a) 
e a N = - 2p (5 + 47V - dy^j 2 , U d 2 > (3 + 57V) 2 , (16b) 

for the symmetric (s) and antisymmetric (a) solutions respectively. The second termination condition Eq. (|15b|) puts 
a constraint on the values the potential parameters Uq and d can take. We shall now give some illustrative examples 
of the first two states N — 1,2. 



A. The N = 1 state 

Let us consider the first symmetric state, N = 1. Here the eigenvalue is e\ = —4^2 (7 — dy/Uo) 2 \ and we must 
ensure Uod 2 > 49 to guarantee Eq. (|15a|) is satisfied. The second termination condition Eq. (|15b|) can be re-written 
as the 2x2 matrix 



p-qi 1 + P 
a p — <?2 + Oi 



where 



<Zi = 0, 



q 2 = 2 + p + 7 



(17) 
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which tells us what pair of values Uq and d can take. For simplicity, in what follows we choose d = 1, and upon 
solving Eq. ([TT|) via root-finding methods we find the special potential strengths 

U = 149.57..., and 595.84... (18) 

Thus, we have found the following wavefunction describes the first symmetric bound state in a hyperbolic double- well 

1>(0. = e*U P/2 + (19) 

We plot Eq. ([T9]) and its associated probability density for the special values Eq. ([15]) in Fig. [2j showing the symmetric 
parity of both wavefunctions as one would expect. We see that the higher the potential strength the tighter the 
confinement and a decrease in nodes from two to zero (although the factor in the wavefunction produces a 
superficial "node" at x = in the second case). Carrying out the same analysis for the first antisymmetric state, we 
find similar behavior but now with a decrease in nodes from three to one, as can be seen in Fig. [3l 




FIG. 2: (Color online) Plots of two symmetric wavefunctions from the N = 1 state, with d — 1 and (a) Uo — 149.57... (blue 
line), (b) Uo — 595.84... (red line). Inset: associated probability density. 
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FIG. 3: (Color online) Plots of two antisymmetric wavefunctions from the N = 1 state, with d — 1 and (a) Uq — 426.23... (blue 
line), (b) Uq = 1092.80... (red line). Inset: associated probability density. 



B. The N — 2 state 



Proceeding as before, but this time for the N = 2 symmetric state, we immediately find from Eq. (fTBj) the eigenvalue 



^2 (ll — d\/Uo) 2 , subject to the stipulation Uod 2 > 121. The second termination condition Eq. (|15b[) is 



4d 
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equivalent to the 3x3 matrix 
/i - qi l + P 

2a fi-q 2 + a 2(2 + p) = 0, where q x = 0, q 2 = 2 + /3 + 7, 4s = 2(3 + /? + 7)- (20) 
a n — qs + a 



Upon solving Eq. ([20]) we find the special values 

U = 279.14..., 860.32..., and 1740.79..., (21) 
again for d = 1. The form of the symmetric wavefunction follows straightforwardly from Eq. (fl3j) 

m = e* € ^ /2 + ^ + (22) 

In Fig. 3J we plot Eq. (|22j) and its associated probability density for the special values Eq. (|2T|) . It can be seen that 
the wavefunctions are symmetric, with a maximum of four nodes for the shallowest state, to two nodes for the middle 
state, to zero nodes for the deepest state. Again, the factor in the wavefunction creates superficial "nodes" at 
x = in the second and third cases. The equivalent antisymmetric solution is displayed in Fig. [5j illustrating the 
anticipated pattern: a drop from five nodes to three nodes to one node as we increase the potential strength. 




FIG. 4: (Color online) Plots of three symmetric wavefunctions from the N — 2 state, with d — 1 and (a) Uo — 279.14... (blue 
line), (b) Uo — 860.32... (red line), (c) Uo — 2539.74... (green line). Inset: associated probability density. 
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FIG. 5: (Color online) Plots of three antisymmetric wavefunctions from the N — 2 state, with d — 1 and (a) Uq — 642.50... 
(blue line), (b) Uq — 1445.59... (red line), (c) Uq — 1740.79... (green line). Inset: associated probability density. 



Higher states N = 2,3,... can be obtained by following the same recipe. For wavefunctions of a state TV, there 
are TV + 1 solutions with a different Uq for a certain d. The symmetric wavefunction confined in the lowest potential 
strength Uq will have 27V nodes, and the node number will decrease by two each time until the deepest symmetric state 
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has no nodes. The antisymmetric wavefunction confined in the lowest potential strength will have the most nodes, 
2N + 1, which decreases by two each time until the deepest state has just one node. As the potential is symmetric, 
we find as expected the parity of the polynomial solutions changes alternately from even to odd as we increase the 
potential strength and hit successively higher eigenvalues. 



IV. CONCLUSION 



We have reported a family of confining hyperbolic potentials which allow one to transform the one-dimensional 
Schrodinger equation to a confluent Heun equation. One case, describing a hyperbolic double-well, can be further 
reduced such that the wavefunctions associated with bound states can be written in terms of Heun polynomials. 
We expect this work to be of general interest due its simplicity and focus on the well-known double-well problem of 
quantum mechanics, but also to be intriguing to those interested in the use of the exotic but increasingly popular 
Heun differential equation in physics. 
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Appendix A: The second termination condition as a tridiagonal determinant 



The second termination condition, Ajv+i(aO = 0> can be represented as the following tridiagonal determinant 14 



li- qi (1 + 0) 
Na fi-q 2 +a 2(2 + 0) 



(N -l)a fi - q 3 + 2a 















/i - QN-i + (N- 2)a (N - 1)(N -1 + 0) 
2a /i — qjy + (N — l)a 

a 









N(N + 0) 
M - Qn+i + Na 



= 0, 



(Al) 



where q n = (n — l)(n + + 7). 
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